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THE STRUCTURE OF TURBULENCE IN FULLY DEVELOPED PIPE FLOW ‘ 

By John Laufeb 


SUMMARY 

Measurements, principally with a hot-wire anemometer , were 
made in fully developed turbulent flow in a 10-inch pipe at 
speeds of approximately 10 and 100 feet per second. Emphasis 
was placed on turbulence and conditions near the wall. The 
results include relevant mean and statistical quantities, such as 
Reynolds stresses, triple correlations, turbulent dissipation, and 
eneryy spectra. It is shown that rates of turbulent-energy 
production, dissipation, and diflusion have sharp maximums 
near the edge of the laminar sublayer and that there exist a 
strong movement of kinetic energy away from this point and an 
equally strong movement of pressure energy toward it. Finally 
it is suggested that, from the standpoint of turbulent structure, 
the field may be divided into three regions: {!) Wall proximity 
where turbulence production, diflusion, and viscous action are 
all of about equal importance; {2) the central region of the pipe 
where energy diflusion plays the predominant role, and ifl) the 
region between (1) and (2) where the local rate of change of 
turbulent-energy production dominates the energy received by 
diflusive action. 

INTRODUCTION 

The one aspect of turbulent shear flow that stands out 
most prominently is the transport of stream properties by 
turbulent motions. The transfer process is fundamental, for 
it not only shapes the mean-flow field through momentum 
transfer but supplies the mechanism by which turbulent 
motions receive energy from the mean-flow field. The well- 
known phenomenological theories were the first attempts to 
give analytical forms for the transfer mechanism by some 
simple physical considerations and thereby succeeded in, 
maldng predictions about the nature of the mean-velocity 
field. Subsequent experiments, however, have clearly dem- 
onstrated the inadequacies of these theories, and in a sys- 
tematic discussion Batchelor (ref. 1) has pointed out the 
inconsistencies and unreal consequences of the assumptions 
involved. Recently, Rotta (ref. 2) and Tchen (ref. 3) pre- 
sented more extensive and deeper analytical treatments of 
the nonisotropic turbulence problem, and, although their 
results show some agreement with the existing experimental 
findings, they cannot escape some arbitrariness in assuming 
the nature of the transfer mechanism. It is the general 
opinion among investigators in this field that more extensive 
experimental work on the turbulence mechanism is necessary 

1 Supersedes NACA TN 2954, “The Structure of Turbulence in Fully Developed Pipe 1 
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before a satisfactory analytical formulation of the problem 
will be possible. 

In the last few years various experimental investigations 
were carried out in different types of shear flows, It became 
apparent that the dynamic and kinematic processes govern- 
ing such flows may be quite different. The difference be- 
tween the so-called free turbulent flows and flows past solid 
boimdaries has recently been emphasized by the phenomenon 
of intermittency, first noted in a jet by Corrsin (ref. 4) and 
later more clearly recognized and studied in detail in the 
wake of a cylinder by Townsend (ref. 5). It is apparent by 
now that in flows like wakes, jets, and those near the free 
surface of the boundary layers the intermittency is present 
and seems to play a very important^ role in the transfer 
mechanism, while in pipe and channel^flows it is completely 
absent. 

Fairly large amounts of experimental information have 
been gathered about flows with free boundaries (refs. 6, 7, 
and 8); especially extensive is the work by Townsend in a 
turbulent wake (ref. 5). The recent work in two-dimen- 
sional channel flow by the present writer (ref. 9) brought out 
some significant features of flows with solid boundaries, but 
the information was far from being complete. One conse- 
quence of this study was the realization that in order to 
obtain a complete picture, say of the turbulent-energy 
balance, a knowledge of flow conditions in close proximity 
to the wall was of utmost importance. 

Since no similar investigation, with emphasis on the turbu- 
lent structure, has ever been carried out in fully developed 
pipe flow, the present investigation was undertaken. A pipe 
also offered the simplicity of axially symmetric mean flow 
and a turbulent field nonhomogeneous in one direction only. 
At the same time it afforded an experimentally convenient 
setup for obtaining a shear flow of large scale. A diameter 
of 10 inches was chosen as sufi&ciently large for the favorable 
application of hot-wires and for bringing the region near the 
waU within practical reach. To cover the desired range, 
two working Reynolds numbers were chosen, these being 
500,000, based on the diameter and an airspeed at the center 
of 100 feet per second, and 50,000 at a speed of 10 feet per 
second. The higher Reynolds number minimized viscous 
effects while the lower one magnified the extent of the 
predominantly viscous layer near the wall. This turned out 
to be a fortunate choice since experimental difficulties and 
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errors encountered in one case were usually absent in the 
other. 

The purpose was to investigate the nature of turbulence 
and its relation to the mean flow, in particular the rates of 
transfer, diffusion, and dissipation of energy. To this end 
measurements were made of the relevant mean and statis- 
tical quantities, including their spatial distributions from 
within the laminar sublayer to the center and energy spectra 
at several radii. Thus the purpose was to reveal details of 
flow structure that could not be found from studies of mean 
flow alone, and yet stop short of a direct attempt to find the 
nature of the transfer mechanism. It is hoped, however, 
that advances toward the latter goal will have been made 
through the present work. 

The present investigation was conducted under the spon- 
sorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The author wishes 
to express his gratitude to Dr. G. B. Schubauer for his 
constant interest and encouragement throughout this work. 
The valuable discussions with Dr. C. M. Tchen and Mr. 
P. S. Klebanoff are also much appreciated. Thanks are due 
to Miss Z. W. Diehl and Mr. K, D, Tidstrom who carried 
out the numerical computations and preparation of the 
report. 

SYMBOLS 


a 

Do, Dr 


^1, ^2 

Foik^ 

F,o{h) 

(QD)r 


k. 

P 

P. 

(PD)o, {PD)r 


pipe radius, 4.86 in. 

dimensionless turbulent-kinetic-energy dif- 
fusion rates; 

aid 


Do^ 


Ur^ r dr 



^~Ur^rdr 

voltage fluctuations across hot-wires 
fraction of turbulent energy associated 
with ki, cmVsec^ 

fraction of turbulent energy v‘‘ associated 
with ki, cm^/sec^ 

fraction of turbulent energy associated 
with ki, cm®/sec^ 

dimensionless gradient diffusion rate of 
timbulent kinetic energy, 

1 d d 

Ur^ r dr ^ dr 2 

one-dimensional wave number in direction 
of mean flow, cm"^ 
mean pressure at any point in pipe 
mean pressure at exit of pipe 
dimensionless turbiflent-pressure-energy dif- 
fusion rates; 

{PD)o = 7 ^ 7 1 ; ^ ^ 


{Pr)o, (Pr)r 


P 

2 

R 

Rr 


P =a—r 
t 

U 

Uo 

. 

'az 


U, V, V) 
u' , v', w' 
T ", W 

Wo, Wr 


w, 

X 

p 

7T 

P 

<P 


dimensionless turbulent-energy production 
^ rates, 

instantaneous value of pressure fluctuations 
dynamic pressure at pipe center 
Reynolds number based on diameter of pipe 
and velocity at pipe center 
correlation coefficient of -it-fluctuations at 
two points displaced in a radial direction 
coordinate in radial direction; r=0 corre- 
sponds to pipe center 


friction distance parameter, r'Urlv 
time 

mean velocity at any point in pipe 
maximum value of mean velocity 


friction velocity; Ur^ = ~v 





total (mean-plus-fluctuating) velocities in x, 
r, and <p directions, respectively 
instantaneous values of velocity fluctua- 
tions in X, r, and <p directions, respectively 
root-mean-square values of velocity fluctua- 
tions in X, r, and <p directions, respectively 
mean velocities in radial and azimuthal 
directions, respectively 
dimensionless turbulent-energy dissipation 
rates; 



where Ui refers to three velocity fluctua- 
tion components u, v, and w and Xj, to 
three coordinates; repeated indices indi- 
cate summation 

dimensionless direct-viscous-dissipation rate, 

Ur\dr) 

dissipation length parameter 
kinematic viscosity of air 
total (mean-plus-fluctuating) pressure 
air density 
azimuthal coordinate 


ANALYTICAL CONSIDERATIONS 


EQUATIONS OF MOTION OP PIPE PLOW 

The continuity equation and Reynolds equation in cylin- 
drical coordinates for incompressible mean flow have the 
following forms : 

.1 brV 1 bW 
bx'r br ' r 


( 1 ) 




and 

dU „i>U Wi,U_ 
dx dr r 
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TLe boundary conditions are at r=0 

^=0 
dU 


1 ^ 
p dx 


(~ 'u^+— ^ ruv-\-~ ^ 

\da; r dr r d^ / 


^-L u 

da; dr r b<p r 


1 ^ 
p dr 


Id , 1 d 

— ^ ry^d — vw 

r or T o<p 




jj dPT_^y dTF , TF dW , FTF 


da; 


dr 


d , 1 d -2 
dr r 0<p 

where 


r b<p T prb<p 

2^+/v2Tir+-^ 


_ F_2 dTF\ 

^<p) 

UW + 


1 dP 


r^ d<p 


-(- 

\da: 

JF\ 

rV 


( 2 ) 


_2_^ , ^ , 1 _d .1 ^ 

dx^~^dr^ ‘ r dr"^P d^^ 

In a fully developed turbulent pipe flow the conditions are 
(a) F=0 and TF=0 

(c) The velocity field is independent of the coordinate x 
Equations (2) therefore become 

l.dP Id _ , fd‘^U I 1 dU\ 
pda; r dr \dr^ r dr J 

1 dP Id ^ I ^ 

p dr r dr r 

^ d — 2Fw5 

0 = — vw 

dr r 

Integrating the last equation and using the boundary condi- 
tion ?w=0 at r—a it follows that ?w=0 for all values of r. 
Thus the Reynolds equations for the turbulent pipe flow 
reduce to 

(3a) 


dx 


r dr ^ 


1 dP Id “21^ 

p dr r dr r 


(3b) 


Differentiating equation (3b) with respect to x one finds 
d2P/drdx=0. Thus dP/dx is independent of r and equa- 
tions (3 a) and (3b) readily integrate to 


dr 


= 0 


and at ? =a 


uv = 0 
dU _ JJ 2 

V -J— = — t/T 

dr 

Also let P=0 at x=0 and r=a. Then from equation (4a), 
A(x)=0, 

ldP_ 2 2 

a ox a 


and, mtegrating, 


^=-2 U,‘x+C{r) 

P d 


(5) 


From equations (4b) and (5) , P(x) = -- Wx. The equations 
finally become 


i dU I r JJ 2 
uv = v -^—-\ — Ut 
dr a 




— W‘ 


J P 2 ^7 2 



p a 


(6a) 

(6b) 


It is interesting to note, as already pointed out by Kampe 
de Feriet in the case of flows between parallel walls (ref. 10), 
that the shearing stress uv and the mean veloc^ occur 
together in one equation and the normal stresses and 
and the mean pressure occur together in the other. From 
the experimental point of view this has the advantage that 
the equations furnish a method of checking the absolute 
accuracy of the measurements of the Reynolds stresses if the 
mean-velocity and pressure distributions are known. 

ENERGY EQUATIONS 

The momentum equations may be written in the form 
, buiU2 . 1 buiUz , U 1 U 2 


dx dr 


d^ 


^iU2 I d^2^ I 1 d'Ua^s [ U2^— 
dx ’ dr ”^r d^ r 


p dx 

Uz^ 1 d-TT , 

—=“57+ 


/_2 '^2 2 dl ^\ 


— — r (^v—v ^-^-i-A.(x) 

(4a) 

buiUz 1 

dw2'**3 1 1 diia^ I 2u2Uz, 

dx 

^ dr r d^ r 

^ r^'“'^'dr+P(x) 

Ja r 

(4b) 




1 dir 


+ 


pr d<p 
/_2 '** 3,2 buz\ 


( 7 ) 
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Multiplying the three equations by Ui, and 1^3, respec- 
tively, one'gets 

1 Ui^U2 2 bx , 

dx dr ' r ‘ r ~~ p^^dx ‘ 

-[i~-(s)'-eT-Ks')'] 

'buiU2^ , d'W2® I 1 dw2®«3 , U2{u2^~2u^ 2 dir , 

r “~p“’dF+ 



/ d^ Y _ / ^ V _ Jl 

\dx/ \dr/ r^\'d<p/ 



2 d'Uj 

r^ d^ 


1 


duiUs^ dUjUa^ ■ 1 d^ 3 ^ Zu 2 Uz^ _ _ 2 % dr 
dx dr ”^r d<p r ~ p r b<p' 



fdUzV 1 fduzY d^3~ 

\dx/ V dr / r^vd^p/ r^*^r^^^dv?_ 


Introducing the velocity and pressure perturbations 

Ui—U-\-u 

U2—V 

Uz = W 

r— P+P 


averaging, and adding the three equations there is obtained 
the energy equation for the turbulent motion 


—,dU , 1 d /' 

uv-^ T-'>''0\ ■ 

dr r dr \ 




)= 


d d . I' / 4 dr 

r dr^ dr 2 ~^r\r^^ r / 


/d^A/duA 
\dxj/ \dx^/ 


(8) 


where in the last term Ui refers to the three velocity fluctua- 
tion components and Xj, to the three coordinates; the repeated 
indices indicate summation. 

The second term on the right-hand side of the equation 
may be neglected compared with the other terms for all 
values of r, since (a) near the center it approaches zero as 
r^O as may be shown from the continuity equation; and (b) 
everywhere else its order of magnitude where Z> is a 

characteristic length of the order of the pipe radius, is 
smaller than, for instancy that of the last term of the 
equation; that is, vu^j\'^:^vu^lD^, since XVP2<1. 

The equation may now be rewritten as 


— dU , I d 

uv T~ 

dr r dr 


I 




■?)- 


V d d 

r dr dr 2~~~ 


'-■ (©©)-» <•> 


The energy equation thus obtained has essentially the 
same form as that given by Von Kdrmdn for the case of a 
turbulent flow between two parallel plates (ref. 11). The 
first term corresponds to the rate of production of turbulent 
energy by the action of the shearing stresses. The second 
term represents the rate of energy change due to transfer of 
both kinetic and potential energies by the radial velocity 
fluctuations and is usually referred to as the diffusion term. 
The third expression may be regarded as a gradient type of 
energy diffusion and is important only very near the wall. 
The fourth term expresses the rate of energy dissipation into 
heat by action of the viscosity. 

EQUIPMENT 

WIND TUNNEL 

The investigation was carried out in the experimental 
setup shown in figure 1. A centrifugal blower having a 
capacity of 12,000 cubic feet per minute and powered by a 
15-horsepower constant-speed motor was placed inside a 
large pressure box. The box was octagonal in cross section 
and was 10 feet high and 16 feet long. The air was dis- 
charged from the blower into the space surrounding the 
blower and then passed through a honeycomb, three screens, 
and a large contraction cone. The honeycomb consisted of 
hexagonal cells 2 inches across and 8 inches long. The’ 
screens had 24 meshes per inch and a wire diameter of 0.0075 
inch. The contraction cone, which was made out of hard- 
wood and was circular in cross section, reduced the air pas- 
sage from 48 inches to 18 inches in diameter. The cross 
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section was then further reduced from 18 inches to 10 inches 
in diameter through a 90° elbow. After the elbow a short 
elastic coupling was used to prevent the transmission of- 
vibration from the pressure box to the pipe.. In order to 
minimize any flow irregularities due to the elbow, another 
screen and an 8-inch-long honeycomb were installed in the 
entrance sections of a 25-foot-long steel pipe. Furthermore, 
in order to accelerate the boundary-layer growth, the pipe 
wall was artificially roughened by gluing floor-sanding paper 
to the surface along a length of 2% feet. With this arrange- 
ment it was found from the measured mean-velocity distribu- 
tion at the end of the steel pipe, that is, after an entrance 
length of about 30 diameters, that the flow was fuUy devel- 
oped turbulent. Following the steel pipe, a 16-foot-long 
seamless brass tube was attached having an inside diameter 
of 9.72 inches. This served as the actual test section. 

The speed of the tunnel was regulated by throttling the 
intake of the blower with adjustable vanes and by venting 
the pressure box. In the early stages of the investigation 
it was realized that for the very low velocity pipe experi- 
ments, where the intake of the blower had to be nearly 
closed, the heat dissipated in the pressure box was sufiicient 
to affect the flow in the test section. This difiBculty was 
overcome by supplying the air to the system by a ventilating 
fan placed on the side of the box as shown in figure 1. 

HOT-WIRE EQUIPMENT 

The basic hot-wire equipment used during the investiga- 
tion is described in detail in reference 12. The hot-wires 
were made of platinum drawn by the Wollaston process. 
Later during the investigation platinum-rhodium wires (90 
percent platinum and 10 percent rhodium) were used and 
were found to be very satisfactory. The wire diameters were 
generally 0.0001 inch; only when the noise problem was 
very critical were finer wires (0.00005 inch) used. For the 
measurements of the longitudinal components of the velocity 
fluctuations the length of the wires ranged from 0.01 to 
0.025 inch. 

Special care was taken in building the X-type of wires for 
the measurements of the cross components v and iv. In 
order to minimize wire-length effects and to be able to work 
very near the wall, the size of the wire holder had to be as 
small as possible. With the use of prongs made of fine 
jeweler’s broaches the dimensions of the holder head were 
cut (Town to approximately 0.015 inch by 0.005 inch, the 
wires having a length about 0.025 inch. 

TRAVERSING MECHANISM | 

The traversing mechanism simply consisted of a microm- 
eter screw on which the hot-wire support was fastened. 
The support could be rotated in a plane perpendicular to 
the air flow so that the hot-wire could be adjusted parallel 
to the wall. 

The zero reading of the traversing mechanism (r'=0) 
was found by placing the hot-wire close to the wall (approxi- 


mately 0.01 inch away) and measuring the distance between 
the wire and its image in the polished wall by an ocular 
micrometer. Since the curvature of the wall was small, the 
space between the wire and the wall could be taken as half 
the observed distance. 

PROCEDURE AND RESULTS 

While fully developed turbulent flow was readily obtained 
throughout the 16-foot test length, tedious and time-, 
consuming adjustments had to be made to remove secondary 
effects and obtain axial symmetry. The task was particularly 
difficult at the low Reynolds numbers where the effect of 
temperature gradient was most felt. This as well as dis- 
turbances from the elbow had to be eliminated. After the 
final adjustments, pressure and velocity surveys made along 
the 16-foot test length showed that the velocity field at each 
section was the same and symmetry about the axis had been 
obtained. All the final measurements were made inside the 
tube 2 to 4 inches from the exit. 

MEASUREMENT OP MEAN VELOCITY AND PRESSURE 

The experiment was conducted at two Reynolds numbers, 
50,000 and 500,000, corresponding to maximum mean veloci- 
ties of approximately 10 and 100 feet per second, these 
velocities being varied slightly in order to compensate for 
the daily changes in air viscosity and density. In order to 
cover the wide range of static and dynamic pressures, an 
inclined manometer 5 feet in length was used. With benzol 
and a slope of 10 to 1 the manometer was sensitive to 0.005 
centimeter of water by direct reading. With a traveling 
microscope and separate scale the sensitivity was increased 
to 0.0001 centimeter of water. 

The pressure distribution in the direction of the flow was 
measured through pressure taps located every 2 feet along 
the brass tube. The results are given in figure 2. 
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Small total-head and static tubes were made of 0.04-inch- 
diameter nickel tube stock with 0.003-inch wall thickness. 
The tip of the total-head^.tube was flattened to an opening of 
0.006 inch. The static tube was placed approximately 0.4 
inch above the total-head tube. This arrangement was used 
for measuring the complete velocity distribution in the low 
Reynolds number case. In the vicinity of the wall a cor- 
rection had to be made because of the large level of velocity 
fluctuations. The correction, having the form 


was of the order of 5 percent or less. 




Figure 4. — Mean-velocity distribution near wall. Dashed lines are 
computed from pressure-drop measurements. 


For the high Reynolds number a hot-wire was used to 
explore the mean-velocity distribution near the wall. Here 
again a correction for velocity fluctuations had to be made 
because of the nonlinear behavior of the hot-wire. This was 
accomplished by an approximate graphical method using 
the known static response curve of the wire (voltage against 
velocity) and the known root-mean-square value of the 
voltage fluctuations. The maximum correction was about 
10 percent, the correct mean velocity being higher than the 
observed. At the low Reynolds number the velocity profile 
determined by the hot-wire was unreasonably low in the 
vicinity of the wall after applying the correction. Except 
for this, the agreement between the hot-wire results and 
those obtained with the total-head tube was good. 

The measurements are presented in figures 3 and 4. In 
figure 4 the dashed lines indicate the wall velocity gradients 
computed from the pressure drop. The agreement with the 
measured values is satisfactory. 



O .1 .2 .3 .4 .5 .6 .7 .8 .9 i.O 

r'/a 


Figure 5. — -u' distribution. 
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Figure 8. — Reynolds shearing stress and double-correiation-coefficient 
distributions. Curves calculated from measured dUldr, u', and v'. 


MEASUREMENTS OF TURBULENCE LEVELS AND SHEARING STRESS 

Tlie three components of the velocity fluctuations u', 
and w' and the turbulent shearing stress uv were obtained by 
standard techniques described in reference 6. 

Figures 5 to 8 give the detailed results. In figure 8 the 
solid lines represent distributions calculated by equation 
(6a) using the independently measured mean-velocity 
g.’adients and u' and v\ It is seen that for the high Reynolds 
number the agreement between these and the directly 
measured points is very good while for the low Reynolds 
number the measured points are somewhat higher. It 
should be mentioned that lower accuracy of all the measure- 
ments in the low Reynolds number flow is to be expected 
mainly because of the difliculty of forming consistent time 
averages because of the inherent low-frequency fluctuations. 
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Figure 10. — Quadruple-correlation coefficients. 


MEASUREMENTS OF TRIPLE AND QUADRUPLE VELOCITY CORRELATIONS 

Measurements of the triple and quadruple velocity corre- 
lations are given in figures 9,10, and 1 1 . The basic technique 
described in reference 12 was essentially adopted for these 
measurements. The schematic diagram of the electronic 
circuit is given below : 
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The two hot-wire signals ei and €2 are fed into two identical 
compensated amplifiers after which ei is squared. The pre- 
liminary squaring circuit uses the nonlinear characteristics 
of two balanced triodes as employed in Townsend’s circuit 
(ref. 13). The output of the mixing circuit, which is a 
simple resistance network (ref. 12), gives simultaneously the 
sum and difference of the inputs. Squaring circuits 1 and 2 
contain a series of diodes properly biased to give the square 
of the inputs. If the signals before the mixing circuit are 
equalized the ratio meter reads directly the correlation 

coefficient 61^62 

It is seen that in order to obtain the triple velocity correla- 
tions by this method the mean fourth power is necessary. 
This could have been avoided, of course, if instead of using 
the ratio meter the outputs of the two squaring circuits had 
been recorded separately and then subtracted. It was felt, 
however, that the method adopted gave more consistent 
results. The mean fourth power was obtained simply by 
feeding the amplified hot-wire signal into a squaring circuit 
and then squaring the output again by reading the fina l 
output on a thermocouple meter. 

Unfortunately, for the lower Reynolds number some in- 
consistency in the measurements was found using the above 
technique. In comparing the values of the double- 
correlation coefficient uvlu'v' obtained by this setup (without 
the preliminary squaring circuit) with those obtained by the 
conventional method, they were found to be from 20 per- 
cent to 30 percent lower everywhere except near the wall. 
At the high Reynolds number the agreement was very 
satisfactory. Although, because of time limitations, it was 



Figure 11 . Triple-velocity-correlation distributions. 


not possible to trace definitely the cause of this discrep- 
ancy, it is believed to be due to a difference in the low- 
frequency phase shift of the two amplifiers. Consequently, 
the large difference between the triple -correlation dis- 
tributions at the two Reynolds numbers in figure 11 should 
be considered to be due to experimental error. 

MEASUREMENTS OF VARIOUS DISSIPATION TERMS 

The expression for the rate of energy dissipation as it 
appears in the turbulent-energy equation (9) has the fprm 




The first three terms were obtained by the differentiation 
method introduced by Townsend (ref. 13). By making the 
assumption 


he writes 


dx U dt 


( 11 ) 




Thus by electronically differentiating the hot-wire signal 


eccu (or v, or w), 



may be easily obtained. 


A recent 


paper by Lin (ref. 14) shows that this assumption is valid 
if no mean-velocity gradient exists and (ulUy<^l. In a 
shear flow he gives an additional condition 


jj bu bU 
bx br 


In the present measurements these conditions were found to 
be satisfied with the exception of a region inside the laminar 
sublayer. Also during the course of a boundary-layer 
investigation at the NBS an experimental verification of the 
validity of this method was made at 0.05 times the boundary- 
layer thickness. 

The fourth and seventh terms of equation (10) were 
obtained essentially by a method first suggested by Taylor 
(ref. 15), in which 


7? - . L 



as r-»0, where u and u(J) denote fluctuations at distance r 
apart. Thus by measuring the correlation coefficient Rr for 


small values of 



can be calculated. 


During the 


present investigation the accuracy of the method was 
greatly improved by adopting the following technique: 
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The above equation may be rewritten as 

Rr 


2u'u'(j^ 2u'u'(^ 


i_ /^Y 

2i? W/ 


Since, for small values of r, u^=v?(J), one has 


<r = 




2u^ 


1_ 2 
W/ 


Several values of <r can be easily measured by placing two 

/duY 

hot-wires at various small distances v apart. Then j 

nay be calculated from the slope of a straight line in a plot 
of <r against P. 

The remaining four terms of equation (10) could also have 
been obtained with the method described above. However, 
the wire arrangement necessary for such a technique made 
its application impractical. Comparing the first three terms 
it was found that they satisfy the isotropic relations fairly 
well except near the wall (figs. 12 and 13). The other two 
measured terms, while not too different from the first three in 

the center region, are considerably higher ^especially Gr)) 

near the wall. This, of course, is not surprising since dis- 
sipation lengths or microscales in the radial direction are 
expected to be smaller because of the presence of the wall. 
For want of a better procedure it was assumed that mean- 
square derivatives with respect to a given direction sepa- 
rately satisfied the isotropic relations 




dr 7 \ dr 7 






1 /d^Y—^ /d-^Y—^ /dt^Y 

r^Vd^’/ P \d<)?/ r^ \d^/ 



Figure 12. — Distributions of dissipation terms. 72—50,000. 



Figure 13. — Distributions of dissipation terms. i2 — 500,000. 


It is not possible to estimate the accuracy of this apump- 
tion at particular points. However, one can determine the 
error made in the total energy dissipation of the pipe cross 
section. Integrating the energy equation (9) across the pipe 
there is obtained 

( 12 ) 

<2^ Jo dr a Jo 

which states that the total energy produced at a given sec- 
tion is dissipated at the same section because of the homo- 
geneity of the field in the ^-direction. Equation (12) was 
found to be satisfied within 10 percent. It is remarked that 
if the isotropic relation between all terms had been assumed 
and the dissipation calculated in the usual way using 


w=n.(gj 

the right-hand side of equation (12) would have been smaller 
by approximately a factor of 2.5. It is still an open ques- 
tion, however, how accurate the approximation is from point 
to point. 

ENERGY-SPECTRUM MEASUREMENTS 

The amplified hot-wire signal was fed into a Hewlett- 
Packard wave analyzer with a frequency range of 10 to 16,000 
cycles per second. The analyzer selectivity characteristics 
were obtained by calibration; two fixed band widths were 
chosen, their effective values being 11 and 42 cycles per 
second approximately. The output of the analyzer was fed 
into a thermocouple circuit. In order to make the measure- 
ments independent of the amplifier frequency response the 
thermocouple output readings of the hot-wire signal were 
matched by a known sine-wave input. 

The measurements are presented in figures 14, 15, and 16. 
Wire-length corrections were applied only to the li'-spectra, 
using the method described in reference 16. 










10 


REPOKT 1174 — ^NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 






13 


THE STRUCTURE OF TURBULENCE 
DISCUSSION 

As already pointed out in the “Introduction,” the present 
understanding of the turbulence phenomenon is not com- 
plete enough to be able to attack directly the problem con- 
cerning the nature of the turbulence transfer mechanism 
that establishes a stable, nondecaying turbulence field such 
as the pipe fiow. The principal aim of the present work 
was an attempt to obtain an over-all picture of the turbulence 
structure without trying to understand the detailed mecha- 
nism responsible for this structure. To this end the energv 
equations of both the mean and turbulent flow serve as a 
useful guide. The former in effect expresses the relative 
magnitude of the mean-flow energy loss to the turbulent 
field as compared with losses due to dissipation by direct 
molecular action. The latter gives a relation between the 
different forms of turbulent-energy rates such as rates of 
production, dissipation, and diffusion. Unfortunately it is 
not possible to obtain from these equations any explicit in- 
formation on how these various energy changes take place. 
However, it is hoped that once their relative importance in 
different regions of the turbulent field is better understood 
one might attack the ultimate problem of the turbulence 
mechanism with greater confidence. 

A different, but already well-established, method of ap- 
proach was also tried. Instead of investigating the different 
rates of changes of the turbulent energy at a given point of 
the field, its spectral distribution in the wave-number space 
was examined. The merits and shortcomings of this method 
will be discussed later. 



Figure 17. Comparison of rate of turbulent-energy production with 
direct-viscous-dissipation rate near wall. 


IN FULLY DEVELOPED PIPE FLOW 

MEAN-ENERGY BALANCE 

Multiplying the integrated momentum equation (6a) by 
the mean-velocity gradient dUjdr, the following equation is 
obtained: 



Thus the energy available because of the pressure drop 
along the pipe is partly converted into turbulent energy and 
is partly directly dissipated. The two terms on the right- 
hand side of the equation are compared in a nondimensional 
form near the wall in figure 17. There are two points that 
should be noted here: (a) The bulk of the direct viscous dis- 
sipation takes place in a very narrow region, r* <^15; (b) the 
position where the laminar shearing stress is equal to the 
turbulent shearing stress (viscous dissipation equal- to turbu- 
lence production) is found to be approximately at the same 
point where the maximum amount of energy is produced 
(/•*=: 11.5). This point is usually referred to as the edge of 
the laminar sublayer. It is seen that not only is the bulk of 
the energy taken from the mean flow directly dissipated but 
a considerable portion of the total turbulence production also 
takes place here. 

It is quite apparent from this picture that, in order to 
obtain the complete picture of the turbulent-energy balance, 
conditions near and within the laminar sublayer have to be 
known. This, of course, is a very difficult task from the 
experimental point of view. 

TURBULENT-ENERGY BALANCE 

From the measurements presented in the previous sections 
all the terms of the turbulent-energy equation (eq, (9)) can 
be calculated with the exception of the pressure-energy 
diffusion. Unfortunately, since this term is small every- 
where except near the wall, its determination from equation 
(9) is very inaccurate. Nevertheless, it is possible to obtain 
in a qualitative way an over-all picture of the complete 
turbulent-energy balance at a cross section. Since the 
experimental problems and errors are different in the two 
Reynolds number flows they will be discussed separately. 

low Reynolds number flow. — In order to be able to study 
the flow conditions very close to the wall it is necessary to 
produce as thick a viscous layer as possible. This may be 
done by carrying out the measurements at very low mean 
speeds. Certain compromises, however, have to be made. 
First it is much more difficult to establish good flow condi- 
tions at low speeds, and second it is more difficult to carry 
out time-averaging processes in the measurements of statis- 
tical quantities. On the other hand the absolute magnitude 
of the dissipation terms may be better established since 
problems such as wire-length effects, amplifier frequency 
response, and noise are not so critical. 

Figures 18 and 19 give the distribution of all the energy 
terms in the low Reynolds number flow. It should be noted 


« 
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Figure 18. — Turbulent-energy balance. = 50,000. Dashed lines 

are estimated. 



; r'/o 

Figure 19. — Turbulent-energy balance near wall. 22 = 50,000. Num- 
bers in parentheses along ordinate correspond to subscript o, as in 
figure 18. Dashed lines are estimated. 

that in forming dimensionless quantities for the coordinates, 
the characteristic length a was chosen in figure 18 and 
vjUr was used in figure 19 representing conditions near the 
wall. The correspondence between respective coordinates 
is indicated parenthetically in figure 19. 

The following interesting conclusions may be drawn: 

(1) Throughout the whole cross section, with the exception 
of the center region, the rate of energy production at a point 


is approximately balanced by the rate of energy dissipation. 

(2) All the various energy rates reach a sharp maximum 
near the edge of the laminar sublayer. 

(3) This edge appears to be also the region from which 
the turbulent kinetic energy is diffusing both toward the 
pipe center and toward the wall and toward which the 
pressure energy is transported. 

There is some question about the direction of the pressure 
diffusion in the center region of the pipe. It was mentioned 
in the previous section that the measured triple velocity 
correlations are believed to be too low in this region. Assum- 
ing that thei^ dimensionless value is approximately the same 
as the values found in the high Reynolds number flow — an 
assumption which holds approximately true for all other 
measured dimensionless statistical quantities — kinetic-energy 
diffusion is estimated as shown with a dashed line in figure 18. 
It is believed that the corresponding pressure-diffusion 
distribution (dashed line) comes closer to the actual picture 
than that obtained from the directly measured points (solid 
line) . This would indicate that the direction of the pressure 
diffusion is toward the laminar sublayer and its value near 
the center is close to zero. 

High Reynolds number flow. — As already pointed out, 
accurate measurements of the various dissipation terms in 
the case of high Reynolds number flow become much more 
difficult, since, because of the extent of the energy spectrum 
to high frequencies, amplifier-tube noise and wire-length 



Figure 20. — Turbulent-energy balance. ^=500,000. Dashed lines 
are estimated. 
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effects become increasingly more critical. Although all pre- 
cautions were made to minimize these effects, the calculated 
rate of energy dissipation is believed to be too small. From 
the energy-spectrum measurements it was possible to infer 
that considerable dissipation takes place at frequencies as 
high as 30,000 cycles per second where the response of the 
compensated amplifier ceases to be linear. From the infor- 
mation gathered in the low Reynolds number flow, an 
attempt was made to estimate the dissipation. The foUow- 
«ing assumptions were made: (1) In the vicinity of the wall 
where similarity with respect to Reynolds number was found 
(this will be discussed in detad later) the dissipation can be 
obtained directly from the low Reynolds number measure- 
ments; (2) the errors in percentage are the same across the 
pipe. This was found to be fairly closely true when the 
dissipation measurements were repeated using various high- 
frequency cut-off filters. With these assumptions and with 
equation (12) the estimated dissipation was obtained. 
Figures 20 and 21 show the directly measured and estimated 
values. The figures also indicate that the picture of the 
energy balance is similar to that obtained in the low Reynolds 


number flow. 
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Figure 21. — Turbulent-energy balance near wall. 22=500,000. 

Dashed lines are estimated. 

ENERGY-SPECTRUM CONSIDERATIONS 

In the study of isotropic turbulence the concept of the 
energy spectrum proved to be a useful and a convenient one 
both from an analytical and a physical point of view. By 
considering the distribution of the turbulent energy in the 
wave-number space it is possible to give in many instances a 
simple picture of some basic physical processes taking place 
in the field. An example is the commonly accepted picture 
of the energy transfer from smaller to larger wave numbers. 
The utility of the spectrum function is much more re- 
stricted in a noniso tropic turbulence field where it cannot be 
contracted into a scalar, as in the isotropic case, but has to 
be treated as a tensor quantity. A restrictive factor in all 
cases, and particularly in nonisotropic turbulence, is the fact 
that measurement is possible of only the so-called Taylor, or 


time, spectrum. This means that, even if Taylor's hypoth- 
esis (eq. (11)) is fulfilled, only a spatial energy distribution 
over a surface Constant can be measured. This, of 
course, prohibits detailed knowledge of the point-by-point 
energy distribution in the wave-number space. Neverthe- 
less, it is possible to draw some conclusions from this type of 
measurements. 

The measured 'li'-spectra (fig. 14) indicate that the spectral 
distributions obtained not too close to the wall manifest 
similar behavior over a wide range of wave numbers. Specifi- 
cally, they vary as the —5/3 power of the wave number k\ 
over a considerable range. This is clearly indicated in 
figure 22 where the function is shown to be a 

constant over the range l<fci<24. The two sets of points 
shown in this figure indicate the magnitude of the length 
corrections applied to the directly measured values. This, 
of course, is the same type of behavior as that of the equilib- 
rium range of the energy spectrum in an isotropic field first 
predicted by Kolmogoroff (ref. 17). One may infer there- 
fore that for sufficiently large turbulent Reynolds number 
flows (ot'X/j^>200 for all measurements) and for the mean- 
velocity gradients not too large there exists a wide range of 
wave numbers in which the energy represented by the 
ii'-component is transferred from smaller to larger wave 
numbers without being significantly influenced by the tur- 
bulent-energy production mechanism or by viscous dissipa- 
tion. This apparently is true, even though local isotropy 
does not exist in this range as will be seen later. It is further 
seen that this range varies with r' fa and at r7a= 0.0082, 
where the mean-velocity gradient is already large, the 
spectrum shows a different distribution. It might be 
mentioned that in this case there is a rather wide wave- 
number range where the spectrum varies closely as as 
predicted by Tchen (ref. 3). 
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Figure 22. — Equilibrium range of u'-spectrum, 22 — 500,000; 
r7a=0.28. 


The conclusions one may draw from the measurements of 
the and w'-spectra are significantly different. In order 
to clarify this difference a comparison is made between the 
©'-spectrum measured at the center of the pipe and that 
calculated from the measured w'-spectrum using isotropic 
relations (fig. 23), It is seen that there is a large energy 
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deficiency in the low wave-number range, while the energy 
content of the higher wave numbers is much larger. Further- 
more, the —5/3-power-law type of distribution is completely 
missing. Similar statements may be made for the y'-spectra 


(and w'-spectra) at points other than the center of the pipe. 
An explanation of such behavior is, of course, difficult 
mainly because the one-dimensional nature of the measure- 
ments conceals the detailed spatial energy distribution. 


Measured 

Calculated 


Figure 23. — Comparison of measured and calculated v' -spectrum, it =600,000; j'7a=1.0 
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However, it is hard to conceive that the geometry of the 
field alone would be responsible for such a behavior. Figure 
23 may suggest that wave numbers receive a large portion 
of their v' energy not only by means of the usual transfer 
mechanism from smaller to larger wave numbers. By con- 
sidering the energy equations for the turbulent components 
separately, the term most likely to influence the -spectrum 
in the wave-number range considered is the pressure term 
transferring energy from one velocity component to the 

other. Its form p ^ suggests that it is associated with wave 

numbers larger than those where the bulk of the turbulent 
energy is produced and smaller than the ones where viscous 
dissipation is important. This term therefore may very well 
be responsible for the larger v' energy present in the wave- 
number range in question. It is difficult to draw any con- 
clusion about the role played by the diffusion terms; in 
general they are thought of as representing a low-frequency 
phenomenon and therefore would influence the w'-spectrum 
in the low wave-number range only. 

It should be emphasized that the above discussion is 
merely speculative in nature and a considerable amount of 
experimental work especially in connection with the pressure 
terms is necessary to be able to give a quantitative picture 
of the energy balance spectrumwise. 

GENERAL CONSIDERATIONS 

On the basis of the presented measurements the flow field 
in the pipe may be divided into three regions exhibiting 
different behavior from the point of view of turbulence 
structure. 

Wall- proximity range. — Measurements near the wall indi- 
cate that the well-known waU-proximity law of Prandtl 
(ref. 18) for the mean velocity may be extended to the 
fluctuating-velocity field also. Thus by using Ur and vjUr 
: as the characteristic velocity and length parameters, the 
, various velocity distributions become independent of the 
Keynolds number in the approximate range 0 <T r* < 30 
(figs. 24 to 26). 

The various energy rates as they appear in the equation 
of turbulent energy are found to be of equal relative im- 
portance; their magnitudes are much larger than those in 
other regions of the turbulent field and they therefore play 
a dominant role in the energy balance over the entire field. 
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Figure 24. — Mean-velocity distribution near wall. 
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Figure 25. — u' , v' , and w' distributions near waU. . 



Figure 26. — Turbulent shearing stress and double-correlation coeffi- 
cient near wall. 


Center portion of pipe. — Since the direct effect of viscosity 
on the turbulent field is negligible in the center portion of 
the pipe, it may be expected that the distributions of the 
fluctuations expressed relative to the characteristic velocity 
Ur are independent of the Reynolds number. Figures 5, 7, 
and 8 show this to be the case. 

This region is^further characterized by the fact that it 
receives a large portion of its turbulent energy by diffusive 
action. 

Intermediate region. — The intermediate range extends 
from r*«100 to r'/a~0.1. By considering the turbulent- 
energy balance in this region (fig. 21) it is seen that the rate 
of energy diffusion is much smaller than that of production 
or dissipation. This means that the energy produced here 
is locally dissipated. It should be noted that this is one of 
the implied assumptions of the mixing-length theories. 
Another assumption, namely, the existence of statistical 
similarity between the velocity components, is, however, 
difficult to accept in view of the discussion in the section 
“Energy-Spectrum Considerations.” Thus, although the 
gradient type of momentum transfer involved in mixing- 
length theories has more experimental support in this region 
than in other portions of the field, its use is not completely 
justifiable. 





18 


REPOBT 1174 — NATIONAL ADVISORT 

SUMMARY OF RESULTS 

Fully developed flow in a large pipe was found to provide 
a very useful medium in which to study the structure of 
turbulence in shear flow. Embodied in the term “structure 
are the interactions between turbulent motions and mean 
flow, and the various transfers and movements of energy 
from point to point and from mean fl’ow through the spectrum 
of turbulent motions. 

The following are the major results: 

1. The importance of a detailed knowledge of conditions 
in the close proximity of the wall was demonstrated. 

2. Using the similarity parameters U-r and vjUr the flow 
field in this region was shown to be independent of the 
Reynolds number. 

3 The various turbulent-energy rates, such as produc- 
tion, diffusion, and dissipation, were found to reach a sharp 
maximum at the edge of the laminar sublayer (r*»12), 
their magnitude being of equal relative importance but much 
larger than those in other regions of the pipe cross section. 

4. It was found that there exist a strong transfer of kinetic 
energy away from the edge of the laminar sublayer and an 
equally strong movement of pressure energy toward it. 

5. In the center region of the pipe, the characteristic 
length and velocity parameters were shown to be a and Ur. 

6. In the region of large mean-velocity gradients but out- 
side of the dissipative region (between r* » 100 and r'ja « 0.1) 
energy diffusion was found to be small compared with the 
turbulent-energy production. 

7. The spectrum measurements indicated that in regions 
where the mean-velocity gradients are not too large the 
ii'-spectra vary as the — 5/3 power of the wave number over 
a considerable wave-number range. 

8. A considerably different distribution of the v'- and w - 
spectra was measured and an explanation of such a behavior 
was attempted. 

National Bureau of Standards, 

Washington, D. C., October 28, 1962. 
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